Implementation of a Mortar Mixed Finite Element Method using a Multiscale Flux Basis

This paper provides a new implementation of a multiscale mortar mixed finite element method for second order elliptic problems. The algorithm uses non-overlapping domain decomposition to reformulate a fine scale problem as a coarse scale mortar interface problem, which is then solved using an iterative method. The original implementation by Arbogast, Pencheva, Wheeler, and Yotov, Multiscale Model. Simul. 2007, required solving one local Dirichlet problem on each subdomain per interface iteration. We alter this implementation by forming a multiscale flux basis. This basis consists of mortar functions representing the individual flux responses for each mortar degree of freedom, on each subdomain independently. The computation of these basis functions requires solving a fixed number of Dirichlet subdomain problems. Taking linear combinations of the multiscale flux basis functions replaces the need to solve any Dirichlet subdomain problems during the interface iteration. This new implementation yields the same solution as the original implementation, and is computationally more efficient in cases where the number of interface iterations is greater than the number of mortar degrees of freedom per subdomain. The gain in computational efficiency increases with the number of subdomains.